Combining LIANA and Tensor-cell2cell to decipher cell-cell communication across multiple samples.

In recent years, data-driven inference of cell-cell communication has helped reveal coordinated biological processes across cell types. Here, we integrate two tools, LIANA and Tensor-cell2cell, which, when combined, can deploy multiple existing methods and resources to enable the robust and flexible identification of cell-cell communication programs across multiple samples. In this work, we show how the integration of our tools facilitates the choice of method to infer cell-cell communication and subsequently perform an unsupervised deconvolution to obtain and summarize biological insights. We explain how to perform the analysis step by step in both Python and R and provide online tutorials with detailed instructions available at https://ccc-protocols.readthedocs.io/. This workflow typically takes ∼1.5 h to complete from installation to downstream visualizations on a graphics processing unit-enabled computer for a dataset of ∼63,000 cells, 10 cell types, and 12 samples.

Defining and benchmarking open problems in single-cell analysis.

With the growing number of single-cell analysis tools, benchmarks are increasingly important to guide analysis and method development. However, a lack of standardisation and extensibility in current benchmarks limits their usability, longevity, and relevance to the community. We present Open Problems, a living, extensible, community-guided benchmarking platform including 10 current single-cell tasks that we envision will raise standards for the selection, evaluation, and development of methods in single-cell analysis.

Network integration of thermal proteome profiling with multi-omics data decodes PARP inhibition.

Complex disease phenotypes often span multiple molecular processes. Functional characterization of these processes can shed light on disease mechanisms and drug effects. Thermal Proteome Profiling (TPP) is a mass-spectrometry (MS) based technique assessing changes in thermal protein stability that can serve as proxies of functional protein changes. These unique insights of TPP can complement those obtained by other omics technologies. Here, we show how TPP can be integrated with phosphoproteomics and transcriptomics in a network-based approach using COSMOS, a multi-omics integration framework, to provide an integrated view of transcription factors, kinases and proteins with altered thermal stability. This allowed us to recover consequences of Poly (ADP-ribose) polymerase (PARP) inhibition in ovarian cancer cells on cell cycle and DNA damage response as well as interferon and hippo signaling. We found that TPP offers a complementary perspective to other omics data modalities, and that its integration allowed us to obtain a more complete molecular overview of PARP inhibition. We anticipate that this strategy can be used to integrate functional proteomics with other omics to study molecular processes.

Multiscale networks in multiple sclerosis.

Complex diseases such as Multiple Sclerosis (MS) cover a wide range of biological scales, from genes and proteins to cells and tissues, up to the full organism. In fact, any phenotype for an organism is dictated by the interplay among these scales. We conducted a multilayer network analysis and deep phenotyping with multi-omics data (genomics, phosphoproteomics and cytomics), brain and retinal imaging, and clinical data, obtained from a multicenter prospective cohort of 328 patients and 90 healthy controls. Multilayer networks were constructed using mutual information for topological analysis, and Boolean simulations were constructed using Pearson correlation to identified paths within and among all layers. The path more commonly found from the Boolean simulations connects protein MK03, with total T cells, the thickness of the retinal nerve fiber layer (RNFL), and the walking speed. This path contains nodes involved in protein phosphorylation, glial cell differentiation, and regulation of stress-activated MAPK cascade, among others. Specific paths identified were subsequently analyzed by flow cytometry at the single-cell level. Combinations of several proteins (GSK3AB, HSBP1 or RS6) and immune cells (Th17, Th1 non-classic, CD8, CD8 Treg, CD56 neg, and B memory) were part of the paths explaining the clinical phenotype. The advantage of the path identified from the Boolean simulations is that it connects information about these known biological pathways with the layers at higher scales (retina damage and disability). Overall, the identified paths provide a means to connect the molecular aspects of MS with the overall phenotype.

Integrating single-cell multi-omics and prior biological knowledge for a functional characterization of the immune system.

The immune system comprises diverse specialized cell types that cooperate to defend the host against a wide range of pathogenic threats. Recent advancements in single-cell and spatial multi-omics technologies provide rich information about the molecular state of immune cells. Here, we review how the integration of single-cell and spatial multi-omics data with prior knowledge-gathered from decades of detailed biochemical studies-allows us to obtain functional insights, focusing on gene regulatory processes and cell-cell interactions. We present diverse applications in immunology and critically assess underlying assumptions and limitations. Finally, we offer a perspective on the ongoing technological and algorithmic developments that promise to get us closer to a systemic mechanistic understanding of the immune system.

Complementing Cell Taxonomies with a Multicellular Analysis of Tissues.

The application of single-cell molecular profiling coupled with spatial technologies has enabled charting of cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single-cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single-cell and spatial datasets to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.

Compartments in medulloblastoma with extensive nodularity are connected through differentiation along the granular precursor lineage.

Medulloblastomas with extensive nodularity are cerebellar tumors characterized by two distinct compartments and variable disease progression. The mechanisms governing the balance between proliferation and differentiation in MBEN remain poorly understood. Here, we employ a multi-modal single cell transcriptome analysis to dissect this process. In the internodular compartment, we identify proliferating cerebellar granular neuronal precursor-like malignant cells, along with stromal, vascular, and immune cells. In contrast, the nodular compartment comprises postmitotic, neuronally differentiated malignant cells. Both compartments are connected through an intermediate cell stage resembling actively migrating CGNPs. Notably, we also discover astrocytic-like malignant cells, found in proximity to migrating and differentiated cells at the transition zone between the two compartments. Our study sheds light on the spatial tissue organization and its link to the developmental trajectory, resulting in a more benign tumor phenotype. This integrative approach holds promise to explore intercompartmental interactions in other cancers with varying histology.

Profiling the heterogeneity of colorectal cancer consensus molecular subtypes using spatial transcriptomics.

The consensus molecular subtypes (CMS) of colorectal cancer (CRC) is the most widely-used gene expression-based classification and has contributed to a better understanding of disease heterogeneity and prognosis. Nevertheless, CMS intratumoral heterogeneity restricts its clinical application, stressing the necessity of further characterizing the composition and architecture of CRC. Here, we used Spatial Transcriptomics (ST) in combination with single-cell RNA sequencing (scRNA-seq) to decipher the spatially resolved cellular and molecular composition of CRC. In addition to mapping the intratumoral heterogeneity of CMS and their microenvironment, we identified cell communication events in the tumor-stroma interface of CMS2 carcinomas. This includes tumor growth-inhibiting as well as -activating signals, such as the potential regulation of the ETV4 transcriptional activity by DCN or the PLAU-PLAUR ligand-receptor interaction. Our study illustrates the potential of ST to resolve CRC molecular heterogeneity and thereby help advance personalized therapy.

A network-based transcriptomic landscape of HepG2 cells uncovering causal gene-cytotoxicity interactions underlying drug-induced liver injury.

Drug-induced liver injury (DILI) remains the main reason for drug development attritions largely due to poor mechanistic understanding. Toxicogenomic to interrogate the mechanism of DILI has been broadly performed. Gene coregulation network-based transcriptome analysis is a bioinformatics approach that potentially contributes to improve mechanistic interpretation of toxicogenomic data. Here we performed an extensive concentration time course response-toxicogenomic study in the HepG2 cell line exposed to 20 DILI compounds, 7 reference compounds for stress response pathways, and 10 agonists for cytokines and growth factor receptors. We performed whole transcriptome targeted RNA sequencing to more than 500 conditions and applied weighted gene coregulated network analysis to the transcriptomics data followed by the identification of gene coregulated networks (modules) that were strongly modulated upon the exposure of DILI compounds. Preservation analysis on the module responses of HepG2 and PHH demonstrated highly preserved adaptive stress response gene coregulated networks. We correlated gene coregulated networks with cell death onset and causal relationships of 67 critical target genes of these modules with the onset of cell death was evaluated using RNA interference screening. We identified GTPBP2, HSPA1B, IRF1, SIRT1, and TSC22D3 as essential modulators of DILI compound-induced cell death. These genes were also induced by DILI compounds in PHH. Altogether, we demonstrate the application of large transcriptome datasets combined with network-based analysis and biological validation to uncover the candidate determinants of DILI.

Metabolic Communication by SGLT2 Inhibition.

BACKGROUND: SGLT2 (sodium-glucose cotransporter 2) inhibitors (SGLT2i) can protect the kidneys and heart, but the underlying mechanism remains poorly understood. METHODS: To gain insights on primary effects of SGLT2i that are not confounded by pathophysiologic processes or are secondary to improvement by SGLT2i, we performed an in-depth proteomics, phosphoproteomics, and metabolomics analysis by integrating signatures from multiple metabolic organs and body fluids after 1 week of SGLT2i treatment of nondiabetic as well as diabetic mice with early and uncomplicated hyperglycemia. RESULTS: Kidneys of nondiabetic mice reacted most strongly to SGLT2i in terms of proteomic reconfiguration, including evidence for less early proximal tubule glucotoxicity and a broad downregulation of the apical uptake transport machinery (including sodium, glucose, urate, purine bases, and amino acids), supported by mouse and human SGLT2 interactome studies. SGLT2i affected heart and liver signaling, but more reactive organs included the white adipose tissue, showing more lipolysis, and, particularly, the gut microbiome, with a lower relative abundance of bacteria taxa capable of fermenting phenylalanine and tryptophan to cardiovascular uremic toxins, resulting in lower plasma levels of these compounds (including p-cresol sulfate). SGLT2i was detectable in murine stool samples and its addition to human stool microbiota fermentation recapitulated some murine microbiome findings, suggesting direct inhibition of fermentation of aromatic amino acids and tryptophan. In mice lacking SGLT2 and in patients with decompensated heart failure or diabetes, the SGLT2i likewise reduced circulating p-cresol sulfate, and p-cresol impaired contractility and rhythm in human induced pluripotent stem cell-derived engineered heart tissue. CONCLUSIONS: SGLT2i reduced microbiome formation of uremic toxins such as p-cresol sulfate and thereby their body exposure and need for renal detoxification, which, combined with direct kidney effects of SGLT2i, including less proximal tubule glucotoxicity and a broad downregulation of apical transporters (including sodium, amino acid, and urate uptake), provides a metabolic foundation for kidney and cardiovascular protection.

Predicting disease severity in multiple sclerosis using multimodal data and machine learning.

BACKGROUND: Multiple sclerosis patients would benefit from machine learning algorithms that integrates clinical, imaging and multimodal biomarkers to define the risk of disease activity. METHODS: We have analysed a prospective multi-centric cohort of 322 MS patients and 98 healthy controls from four MS centres, collecting disability scales at baseline and 2 years later. Imaging data included brain MRI and optical coherence tomography, and omics included genotyping, cytomics and phosphoproteomic data from peripheral blood mononuclear cells. Predictors of clinical outcomes were searched using Random Forest algorithms. Assessment of the algorithm performance was conducted in an independent prospective cohort of 271 MS patients from a single centre. RESULTS: We found algorithms for predicting confirmed disability accumulation for the different scales, no evidence of disease activity (NEDA), onset of immunotherapy and the escalation from low- to high-efficacy therapy with intermediate to high-accuracy. This accuracy was achieved for most of the predictors using clinical data alone or in combination with imaging data. Still, in some cases, the addition of omics data slightly increased algorithm performance. Accuracies were comparable in both cohorts. CONCLUSION: Combining clinical, imaging and omics data with machine learning helps identify MS patients at risk of disability worsening.

Collaborative effect of Csnk1a1 haploinsufficiency and mutant p53 in Myc induction can promote leukemic transformation.

ABSTRACT: It is still not fully understood how genetic haploinsufficiency in del(5q) myelodysplastic syndrome (MDS) contributes to malignant transformation of hematopoietic stem cells. We asked how compound haploinsufficiency for Csnk1a1 and Egr1 in the common deleted region on chromosome 5 affects hematopoietic stem cells. Additionally, Trp53 was disrupted as the most frequently comutated gene in del(5q) MDS using CRISPR/Cas9 editing in hematopoietic progenitors of wild-type (WT), Csnk1a1-/+, Egr1-/+, Csnk1a1/Egr1-/+ mice. A transplantable acute leukemia only developed in the Csnk1a1-/+Trp53-edited recipient. Isolated blasts were indefinitely cultured ex vivo and gave rise to leukemia after transplantation, providing a tool to study disease mechanisms or perform drug screenings. In a small-scale drug screening, the collaborative effect of Csnk1a1 haploinsufficiency and Trp53 sensitized blasts to the CSNK1 inhibitor A51 relative to WT or Csnk1a1 haploinsufficient cells. In vivo, A51 treatment significantly reduced blast counts in Csnk1a1 haploinsufficient/Trp53 acute leukemias and restored hematopoiesis in the bone marrow. Transcriptomics on blasts and their normal counterparts showed that the derived leukemia was driven by MAPK and Myc upregulation downstream of Csnk1a1 haploinsufficiency cooperating with a downregulated p53 axis. A collaborative effect of Csnk1a1 haploinsufficiency and p53 loss on MAPK and Myc upregulation was confirmed on the protein level. Downregulation of Myc protein expression correlated with efficient elimination of blasts in A51 treatment. The "Myc signature" closely resembled the transcriptional profile of patients with del(5q) MDS with TP53 mutation.

Spatial transcriptomics of B cell and T cell receptors reveals lymphocyte clonal dynamics.

The spatial distribution of lymphocyte clones within tissues is critical to their development, selection, and expansion. We have developed spatial transcriptomics of variable, diversity, and joining (VDJ) sequences (Spatial VDJ), a method that maps B cell and T cell receptor sequences in human tissue sections. Spatial VDJ captures lymphocyte clones that match canonical B and T cell distributions and amplifies clonal sequences confirmed by orthogonal methods. We found spatial congruency between paired receptor chains, developed a computational framework to predict receptor pairs, and linked the expansion of distinct B cell clones to different tumor-associated gene expression programs. Spatial VDJ delineates B cell clonal diversity and lineage trajectories within their anatomical niche. Thus, Spatial VDJ captures lymphocyte spatial clonal architecture across tissues, providing a platform to harness clonal sequences for therapy.

Mapping cardiac remodeling in chronic kidney disease.

Patients with advanced chronic kidney disease (CKD) mostly die from sudden cardiac death and recurrent heart failure. The mechanisms of cardiac remodeling are largely unclear. To dissect molecular and cellular mechanisms of cardiac remodeling in CKD in an unbiased fashion, we performed left ventricular single-nuclear RNA sequencing in two mouse models of CKD. Our data showed a hypertrophic response trajectory of cardiomyocytes with stress signaling and metabolic changes driven by soluble uremia-related factors. We mapped fibroblast to myofibroblast differentiation in this process and identified notable changes in the cardiac vasculature, suggesting inflammation and dysfunction. An integrated analysis of cardiac cellular responses to uremic toxins pointed toward endothelin-1 and methylglyoxal being involved in capillary dysfunction and TNFα driving cardiomyocyte hypertrophy in CKD, which was validated in vitro and in vivo. TNFα inhibition in vivo ameliorated the cardiac phenotype in CKD. Thus, interventional approaches directed against uremic toxins, such as TNFα, hold promise to ameliorate cardiac remodeling in CKD.

Multicellular factor analysis of single-cell data for a tissue-centric understanding of disease.

Biomedical single-cell atlases describe disease at the cellular level. However, analysis of this data commonly focuses on cell-type-centric pairwise cross-condition comparisons, disregarding the multicellular nature of disease processes. Here, we propose multicellular factor analysis for the unsupervised analysis of samples from cross-condition single-cell atlases and the identification of multicellular programs associated with disease. Our strategy, which repurposes group factor analysis as implemented in multi-omics factor analysis, incorporates the variation of patient samples across cell-types or other tissue-centric features, such as cell compositions or spatial relationships, and enables the joint analysis of multiple patient cohorts, facilitating the integration of atlases. We applied our framework to a collection of acute and chronic human heart failure atlases and described multicellular processes of cardiac remodeling, independent to cellular compositions and their local organization, that were conserved in independent spatial and bulk transcriptomics datasets. In sum, our framework serves as an exploratory tool for unsupervised analysis of cross-condition single-cell atlases and allows for the integration of the measurements of patient cohorts across distinct data modalities.

Microbiome-based risk prediction in incident heart failure: a community challenge.

Heart failure (HF) is a major public health problem. Early identification of at-risk individuals could allow for interventions that reduce morbidity or mortality. The community-based FINRISK Microbiome DREAM challenge (synapse.org/finrisk) evaluated the use of machine learning approaches on shotgun metagenomics data obtained from fecal samples to predict incident HF risk over 15 years in a population cohort of 7231 Finnish adults (FINRISK 2002, n=559 incident HF cases). Challenge participants used synthetic data for model training and testing. Final models submitted by seven teams were evaluated in the real data. The two highest-scoring models were both based on Cox regression but used different feature selection approaches. We aggregated their predictions to create an ensemble model. Additionally, we refined the models after the DREAM challenge by eliminating phylum information. Models were also evaluated at intermediate timepoints and they predicted 10-year incident HF more accurately than models for 5- or 15-year incidence. We found that bacterial species, especially those linked to inflammation, are predictive of incident HF. This highlights the role of the gut microbiome as a potential driver of inflammation in HF pathophysiology. Our results provide insights into potential modeling strategies of microbiome data in prospective cohort studies. Overall, this study provides evidence that incorporating microbiome information into incident risk models can provide important biological insights into the pathogenesis of HF.

Enablers and challenges of spatial omics, a melting pot of technologies.

Spatial omics has emerged as a rapidly growing and fruitful field with hundreds of publications presenting novel methods for obtaining spatially resolved information for any omics data type on spatial scales ranging from subcellular to organismal. From a technology development perspective, spatial omics is a highly interdisciplinary field that integrates imaging and omics, spatial and molecular analyses, sequencing and mass spectrometry, and image analysis and bioinformatics. The emergence of this field has not only opened a window into spatial biology, but also created multiple novel opportunities, questions, and challenges for method developers. Here, we provide the perspective of technology developers on what makes the spatial omics field unique. After providing a brief overview of the state of the art, we discuss technological enablers and challenges and present our vision about the future applications and impact of this melting pot.

Expanding the coverage of regulons from high-confidence prior knowledge for accurate estimation of transcription factor activities.

Gene regulation plays a critical role in the cellular processes that underlie human health and disease. The regulatory relationship between transcription factors (TFs), key regulators of gene expression, and their target genes, the so called TF regulons, can be coupled with computational algorithms to estimate the activity of TFs. However, to interpret these findings accurately, regulons of high reliability and coverage are needed. In this study, we present and evaluate a collection of regulons created using the CollecTRI meta-resource containing signed TF-gene interactions for 1186 TFs. In this context, we introduce a workflow to integrate information from multiple resources and assign the sign of regulation to TF-gene interactions that could be applied to other comprehensive knowledge bases. We find that the signed CollecTRI-derived regulons outperform other public collections of regulatory interactions in accurately inferring changes in TF activities in perturbation experiments. Furthermore, we showcase the value of the regulons by examining TF activity profiles in three different cancer types and exploring TF activities at the level of single-cells. Overall, the CollecTRI-derived TF regulons enable the accurate and comprehensive estimation of TF activities and thereby help to interpret transcriptomics data.

A network medicine approach to study comorbidities in heart failure with preserved ejection fraction.

BACKGROUND: Comorbidities are expected to impact the pathophysiology of heart failure (HF) with preserved ejection fraction (HFpEF). However, comorbidity profiles are usually reduced to a few comorbid disorders. Systems medicine approaches can model phenome-wide comorbidity profiles to improve our understanding of HFpEF and infer associated genetic profiles. METHODS: We retrospectively explored 569 comorbidities in 29,047 HF patients, including 8062 HFpEF and 6585 HF with reduced ejection fraction (HFrEF) patients from a German university hospital. We assessed differences in comorbidity profiles between HF subtypes via multiple correspondence analysis. Then, we used machine learning classifiers to identify distinctive comorbidity profiles of HFpEF and HFrEF patients. Moreover, we built a comorbidity network (HFnet) to identify the main disease clusters that summarized the phenome-wide comorbidity. Lastly, we predicted novel gene candidates for HFpEF by linking the HFnet to a multilayer gene network, integrating multiple databases. To corroborate HFpEF candidate genes, we collected transcriptomic data in a murine HFpEF model. We compared predicted genes with the murine disease signature as well as with the literature. RESULTS: We found a high degree of variance between the comorbidity profiles of HFpEF and HFrEF, while each was more similar to HFmrEF. The comorbidities present in HFpEF patients were more diverse than those in HFrEF and included neoplastic, osteologic and rheumatoid disorders. Disease communities in the HFnet captured important comorbidity concepts of HF patients which could be assigned to HF subtypes, age groups, and sex. Based on the HFpEF comorbidity profile, we predicted and recovered gene candidates, including genes involved in fibrosis (COL3A1, LOX, SMAD9, PTHL), hypertrophy (GATA5, MYH7), oxidative stress (NOS1, GSST1, XDH), and endoplasmic reticulum stress (ATF6). Finally, predicted genes were significantly overrepresented in the murine transcriptomic disease signature providing additional plausibility for their relevance. CONCLUSIONS: We applied systems medicine concepts to analyze comorbidity profiles in a HF patient cohort. We were able to identify disease clusters that helped to characterize HF patients. We derived a distinct comorbidity profile for HFpEF, which was leveraged to suggest novel candidate genes via network propagation. The identification of distinctive comorbidity profiles and candidate genes from routine clinical data provides insights that may be leveraged to improve diagnosis and identify treatment targets for HFpEF patients.